Black-hole horizons as probes of black- hole dynamics II: geometrical insights 
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In a companion paper fjj], we have presented a cross-correlation approach to near-horizon physics in which 
bulk dynamics is probed through the correlation of quantities defined at inner and outer spacetime hypersurfaces 
acting as test screens. More specifically, dynamical horizons provide appropriate inner screens in a 3 + 1 setting 
and, in this context, we have shown that an effective-curvature vector measured at the common horizon produced 
in a head-on collision merger can be correlated with the flux of linear Bondi-momentum at null infinity. In this 
paper we provide a more sound geometric basis to this picture. First, we show that a rigidity property of 
dynamical horizons, namely foliation uniqueness, leads to a preferred class of null tetrads and Weyl scalars on 
these hypersurfaces. Second, we identify a heuristic horizon newslike function, depending only on the geometry 
of spatial sections of the horizon. Fluxes constructed from this function offer refined geometric quantities to be 
correlated with Bondi fluxes at infinity, as well as a contact with the discussion of quasilocal 4-momentum on 
dynamical horizons. Third, we highlight the importance of tracking the internal horizon dual to the apparent 
horizon in spatial 3-slices when integrating fluxes along the horizon. Finally, we discuss the link between the 
dissipation of the nonstationary part of the horizon's geometry with the viscous-fluid analogy for black holes, 
introducing a geometric prescription for a "slowness parameter" in black-hole recoil dynamics. 

PACS numbers: 04.30.Db, 04.25.dg, 04.70.Bw, 97.60.Lf 



I. INTRODUCTION 

In Ref. 01 (paper I hereafter) a cross-correlation method- 
ology for studying near-horizon strong-field physics was out- 
lined. Spacetime dynamics was probed through the cross- 
correlation of timeseries hi nu and h out defined as geometric 
quantities on inner and outer hypersurfaces, respectively. The 
latter are understood as test screens whose geometries respond 
to the bulk dynamics, so that the (global) functional structure 
of the constructed cross-correlations encodes some of the fea- 
tures of the bulk geometry. This is in the spirit of reconstruct- 
ing spacetime dynamics in an inverse-scattering picture. In 
the context of asymptotically flat black-hole (BH) spacetimes, 
the BH event horizon $ and future null infinity provide 
natural test hypersurfaces from a global perspective. However, 
when a 3 + 1 approach is adopted for the numerical construc- 
tion of the spacetime, dynamical trapping horizons H provide 
more appropriate hypersurfaces to act as inner test screen^ 
In the application of this correlation strategy to the study of 
BH post-merger recoil dynamics, an effective-curvature vec- 
tor Kf{v) was constructed QJJ] on % as the quantity h mrL to 
be cross-correlated with h out , where the latter is the flux of 
Bondi linear momentum (dPf /du) (u) at J? + (here, u and v 
denote, respectively, advanced and retarded times^l). In this 
paper we explore some geometric structures underlying and 
extending the heuristic construction in yl of this effective lo- 
cal probe into BH recoil dynamics. 



1 In paper I future outer trapping horizons were denoted by H + to dis- 
tinguish them from past outer trapping horizons W~ occurring in the 
Robinson-Trautman model, extending the study in Q]. 
Cross-correlation of quantities at H and requires the choice of a 

gauge mapping between the advanced and retarded times u and v. This 
time-stretching issue is discussed in paper I. 



The adaptation of geometric structures and tools from 
to BH horizons is at the basis of important geometric devel- 
opments in BH studies, notably the quasilocal frameworks 
of isolated and dynamical trapping horizons iQ-Ht] (see also 
Refs. @> 01). In this spirit, the construction of Kf s (v) on 
the horizon H partially mimics the functional structure of 
the flux of Bondi linear momentum at J^ + . In particular, 
(dPP /du)(u) can be expressed in terms of (the dipolar part 
of) the square of the news function J\f on sections of J? + , 
whereas the definition of Kf (v) involves the (dipolar part of 
the) square of a function J\f constructed from the Ricci scalar 
2 R on sections of H. However, the functions J\f and JV differ 
in their spin-weight and, more importantly, they show a differ- 
ent behavior in time: whereas jV(u) is an object well-defined 
in terms of geometric quantities on time sections S u C J^ + , 
nothing guarantees this local- in-time character of JV(v) [see 
Eq. (21 below]. The latter is a crucial characteristic of the 
news function, so that JV (v) cannot be considered as a valid 
newslike function on H. 

These structural differences suggest that, in spite of the suc- 
cess of Kf in capturing effectively (at the horizon) some 
qualitative aspects of the flux of Bondi linear momentum (at 
null infinity), a deeper geometric insight into the dynamics 
of H can provide hints for a refined correlation treatment. In 
this context, the specific goals in this paper are: i) to justify 
the role of Kf as an effective quantity to be correlated to 
(dPP /du), suggesting candidates offering a refined version; 
ii) to explore the introduction of a valid newslike function on 
H, only depending on the geometry of sections S v C H; iii) to 
establish a link between the cross-correlation approach in |Q]] 
and other approaches to the study of the BH recoil based on 
quasilocal momentum. 

The paper is organized as follows. Section|II]introduces the 
basic elements on the inner screen H geometry and revisits the 
effective-curvature vector of paper I. Aiming at understand- 



ing the dynamics of the latter, a geometric system governing 
the evolution of the intrinsic curvature along the horizon % is 
discussed, making apparent the key driving role of the Weyl 
tensor. In Sec. [Ill] some fundamental results on dynamical 
horizons are discussed, in particular a rigidity structure en- 
abling a preferred choice of null tetrad on T~L. Proper con- 
tractions of the latter with the Weyl tensor lead in Sec. UVIto 
newslike functions and associated Bondi- 1 ike fluxes on H pro- 
viding refined quantities on the horizon to be correlated with 
Bondi fluxes at J? + , as well as making contact with quasilo- 
cal approaches to BH linear momentum. In Sec. [V] our ge- 
ometric discussion is related to the viscous-fluid analogy of 
BH horizons, providing in particular a geometric prescription 
for the slowness parameter P in 0]. Conclusions are pre- 
sented in Sec.[VT] Finally a first appendix gathers the geomet- 
ric notions in the text, whereas a second appendix emphasizes 
the physical relevance of internal horizons when computing 
fluxes along H. We use a spacetime signature (—,+,+,+), 
with abstract index notation (first letters, a, b, c..., in Latin al- 
phabet) and Latin midalphabet indices, k..., for spacelike 
vectors. We also employ the standard convention for the sum- 
mation over repeated indices. All the quantities are expressed 
in a system of units in which c = G = 1. 



H. GEOMETRIC EVOLUTION SYSTEM ON THE 
HORIZON: THE ROLE OF THE WEYL TENSOR 

A. The inner screen H 

Let us consider a BH spacetime (.M, g a b), with associated 
Levi-Civita connection V a , endowed with a 3 + 1 spacelike 
foliation {St}. Let us consider an inner hypersurface H, to be 
later identified with the BH horizon, such that the intersection 
of the slices S t with the world-tube H defines the foliation of 
H by closed spacelike surfaces {St}- We consider an evolu- 
tion vector h a along H, characterized as that vector tangent 
to T-L and normal to the slices {St} that transports the slice St 
onto the slice St+st- The normal plane at each point of St 
can be spanned in terms of the outgoing null vector l a and the 
ingoing vector k a , chosen to satisfy £ a k a = —1. Directions 
of £ a and k a are fixed, though a rescaling freedom remains 
(see FigUJ). In particular, and without loss of generality in our 
context, we can write 121 



h a = t - Ck a 1 



(1) 



so that h a h a = 2C. Therefore: h a is, respectively, spacelike 
if C > 0, null if C = 0, and timelike if C < 0. 

Regarding the intrinsic geometry on St, the induced metric 
is denoted by q a b, its Levi-Civita connection by 2 D a and the 
corresponding Ricci curvature scalar by 2 R. The area form is 
\ = ^fqdx 1 A dx 2 and we will denote the area measure as 
dA = ^Jqd 2 x. The infinitesimal evolution of the intrinsic ge- 
ometry along H, i.e., the evolution of the induced geometry 
q a b along h a , defines the deformation tensor Q^/ [cf. Equa- 
tion (IA3b in AppendixlAl 




FIG. 1. Worldtube H foliated by closed spacelike surfaces {St} as 
the result of a 3 + 1 spacelike foliation {S t }. The evolution vec- 
tor h a (tangent to T-L and normal to {St}) transports the slice St to 
St+st- The normal plane at each point of St can be spanned by the 
outgoing and ingoing null normal vectors £ a and k a or by n a , the 
unit timelike normal to Et, and s a , the spacelike outgoing normal to 
St and tangent to E t (cf. AppendixlAl . 

< ^ > q ab , referred to as the expansion 
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(2) 



where the trace 6^ 
along h a , measures the infinitesimal evolution of the St ele 
ment of area along "H, whereas the traceless shear con- 
trols the deformations of the induced metric (see Eq. ( 1A6I ) in 
AppendixlAl). Here Sh can be identified with the projection on 
T-L of the Lie derivative Ch [see Eq. dA2t and the remark af- 
ter Eq. JA9H , Before reviewing the effective-curvature vector 
Kf , let us discuss the time parametrization of H. 

We recall that jumps of apparent horizons (AHs) are generic 
in 3 + 1 evolutions of BH spacetimes. The dynamical trap- 
ping horizon framework offers a spacetime insight into this 
behavior by understanding the jumps as corresponding to 
marginally trapped sections of a (single) hypersurface bending 
in spacetime, but multiply foliated by spatial hypersurfaces in 
the 3 + 1 foliation {S t } floUTl. In the particular case of 
binary BH mergers this picture predicts, after the moment of 
its first appearance, the splitting of the common AH into two 
horizons: a growing external common horizon and a shrink- 
ing internal common horizon jT^flill . It is standard to track 
the evolution of the external common horizon, the proper AH, 
but to regard the internal common horizon as physically irrel- 
evant. In Appendix IB1 we stress however the relevance of the 
internal horizon in the context of the calculation of physical 
fluxes into the black-hole singularity. 

In Fig. |2] we illustrate this picture in a simplified (spher- 
ically symmetric) collapse scenario that retains the relevant 
features of the discussion. On one side, the relevant outer 
screen boundary (namely, null infinity J^ + ) is parametrized 
by the retarded time u, something explicitly employed in the 
expression of the flux of Bondi momentum in Eqs. (33) and 
(34) of paper I. On the other side, from the 3 + 1 perspective, 
the moment t c of first appearance of the (common) horizon 
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FIG. 2. Carter-Penrose diagram (corresponding, for simplicity, 
to a generic spherically symmetric collapse) illustrating the time 
parametrization of the outer and inner screens. The outer bound- 
ary given by is properly parametrized by the retarded time u, 
whereas an advanced time v runs along inner boundaries, in particu- 
lar the dynamical horizon Ti. Given a 3+1 foliation {S t }, t c denotes 
the time t at which the horizon first appears. For t > t c , £< slices 
intersect multiply the hypersurface Ti, giving rise to internal and ex- 
ternal horizons. On the contrary, the advanced coordinate v provides 
a good parametrization of Ti from an initial v > Vo. 



corresponds to the coordinate time t at which the 3 + 1 fo- 
liation {S t } firstly intersects the dynamical horizon Ti. For 
t > t c , St slices intersect twice (multiply, in the generic case) 
the hypersurface Ti giving rise to the external and internal 
common horizons (cf. Ti in Fig. [2]). Therefore, the time func- 
tion t is not a good parameter for the whole dynamical horizon 
Ti. An appropriate parametrization of this hypersurface T-L is 
given in terms of an advanced time, such as v, parametrizing 
past null infinity J^~. More precisely, (for a spacelike world- 
tube portion of Ti) we can label sections of Ti by an advanced 
time v starting from an initial value Vq corresponding to the 
first v = const null hypersurface hitting the spacetime singu- 
larity, i.e., Ti = \J v>Vo S v . 



B. Effective-curvature vector 

In paper I the effective-curvature vector Kf was intro- 
duced using the parametrization of Ti by the time function 
t associated with the spacetime 3 + 1 slicing. In particular, 
Kf s (t) was defined only on the external part of the horizon 
Ti, for t > t c . We can now extend the definition of Kf to 
the whole horizon Ti (more precisely, to a spacelike world- 
tube portion of it) by making use of its parametrization by the 
advanced time v adapted to the 3+1 slicing of Ti. Given 
a section S v C Ti, we consider a vector transverse to it 



(i.e., generically not tangent to S v ) and tangent to the 3-slice 
E t that intersects Ti at S v (i.e., S v = Ti l~l E t ). Then, the 
component K eS [£](v) is expressed a^| 

K cfi [^}(v) cx -£ (Csi) (AfivtfdA, (3) 

where Sj is the spacelike normal to S v and tangent to £ t , and 

A/»= f 2 R(v')dv' + Af V0 , (4) 

where 2 R is the Ricci curvature scalar on (S v ,q a b) and Af Va is 
an initial function to be fixed. As commented above, in spite 
of the formal similarity with the news function Af(u) at ^ + 
[cf. Equation (34) in paper I), definition (01 does not guarantee 
the local-in-time character of Af(v) since it is expressed in 
terms of a time integral on the past history. 

In order to study the dynamics of Kf s , we consider the 
evolution of the Ricci scalar curvature 2 R along the world-tube 
Ti . In terms of the elements introduced above, the evolution 
of the Ricci scalar curvature 2 R along h a has the form 

S h 2 R - -8^ 2 R + 2 2 D a2 D h o { ^ - 2 A0W , (5) 

where 2 A = q ab2 D a 2 Di, denotes the Laplacian on St- Expres- 
sion (|5j is a fundamental one in our work and it applies to any 
hypersurface Ti foliated by closed surfaces St- Contact with 
BHs is made when Ti is taken as the spacetime event hori- 
zon or as the dynamical horizon associated with the foliation 

C. Geometry evolution on BH horizons 

We briefly recall the notions of BH horizon relevant here 
and refer to Appendix [A] for a systematic presentation of the 
notation. First, the event horizon (EH) S is the boundary 
of the spacetime region from which no signal can be sent to 
J^ + , i.e., the region in M. not contained in the causal past 
J~(J^ + ) of J" + . The EH is a null hypersurface, character- 
ized as $ = dJ~(J? + ) n Ai. Second, a dynamical hori- 
zon (DH) or (dynamical) future outer trapping horizon Ti is a 
quasilocal model for the BH horizon based on the notion of a 
world-tube of AHs. More specifically, a future outer trapping 
horizon Ti is a hypersurface that can be foliated by marginally 
(outer) trapped surfaces St, i.e., Ti = UteK 1 ^ 4 w ^ tn out g°" 
ing expansion = on St, satisfying: i) a future con- 
dition 0w < 0, and ii) an outer condition 5k6™' < 0. In 
the dynamical regime, i.e., when matter and/or radiation cross 
the horizon (namely when 5^6^ ^ 0), the outer condition is 
equivalent to the condition that Ti is spacelike 1 150 There- 
fore, for dynamical trapping horizons we have C > in Eq. 
(Q3 [cf. discussion after Eq. (IA12H . 



3 For avoiding the introduction of lapse functions related to different 
parametrizations of T-L, we postpone the fixing of the coefficient to Sec. lIVI 
We note that a global constant factor is irrelevant for cross-correlations. 

4 This property actually substitutes the outer condition in the DH characteri- 
zation Bl flal of quasilocal horizons. 
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For both EHs and DHs, an important area theorem holds: 
8 h A= <£ c O^dA > 0. In the case of an EH, Hawking's area 
theorem ¥l7l [l8ll guarantees the growth of the area, whereas 
in the case of a DH, the positivity of 5 h A = — <f s CO^dA 
[cf. Equation IA1 II is guaranteed by its spacelike character 
(C > 0) together with the future condition 6^> < 0. 

We make now contact with Eq. (0 and interpret the ele- 
ments that determine the dynamics of 2 R. The growth of the 
area of a BH horizon guarantees the (average) positivity of 
0( h \ This offers a qualitative understanding of the dynami- 
cal decay of 2 R: the first term in the right-hand side drives an 
exponential-like decay of the Ricci scalar curvature. More 
precisely, nonequilibrium deformations of the Ricci scalar 
curvature 2 R in BH horizons decay exponentially as long as 
the horizon grows in area. Regarding the elliptic operators 
acting on the shear and the expansion [second and third terms 
in the right hand side of Eq. (|5]l] they provide dissipative terms 
smoothing the evolution of 2 R. Indeed, in Sec. [V] we will 
review a viscosity interpretation of 6^ and cr^' , in partic- 
ular associating with them respective decay and oscillation 
timescales of the horizon geometry. 

1. Complete evolution system driving 2 R 

A further understanding of Eq. (0 requires a control of the 
dynamics of the shear <j$ , of the expansion 9^ and of the 
induced metric q ab , the latter controlling the elliptic operators 
2 D a2 D b and 2 A. Therefore, we need evolution equations de- 
termining 5 h q ab , 5 h 6 w and 5 h a^: 

i) ShQab- definition of the deformation tensor. The evolution 
of q ab is dictated by aty an d t c f- Equation (0]. 

ii) ShO^ : focusing or Raychadhuri-like equation. The evo- 
lution of 9^ involves the Ricci tensor R ab , i.e., the "trace 
part" of the spacetime Riemann tensor R a bc d, thus introduc- 
ing the stress-energy tensor T ab through Einstein equations. 

iii) ShO- 1 ^: tidal equation. The evolution of is driven 
by the Weyl tensor C a bc d, i-e-, the traceless part of the space- 
time Riemann tensor, thus involving dynamical gravitational 
degrees of freedom but not directly the Einstein equations. 

The structural feature that we want to underline about these 
equations is shared by evolution systems on EHs and DHs, 
although the explicit form of the equations differ in both cases. 
More specifically, whereas for EHs the evolution equations 
for 2 R, q ab , 9^ and form a "closed" evolution system, 
in the DH case additional geometric objects (requiring further 
evolution equations) are brought about through the evolution 

equations Shq a b, 8h,9^ and 8ho-^ . Moreover, an explicit 
dependence on the function C, related to the choice of 3 + 1 
slicing as discussed later [cf. Equation (1131) 1. is involved in the 
DH case. For these reasons, and for simplicity, in the rest of 
this subsection we restrict our discussion to the case of an EH, 
indicating that the main qualitative conclusion also holds for 
DHs, whose details will be addressed elsewhere. 

The EH $ is a null hypersurface generated by the evolution 
vector h a , a null vector in this case: h a = i a . The null gener- 



ator £ a satisfies a pregeodesic equation £ C V c £ a = K^'l a [see 
Eq. (1A9I > for the expression of the nonaffinity parameter k"]. 
Choosing an affine reparametrization such that l a is geodesic, 
i.e., «W = 0, the evolution equations for 2 R, q ab , and 
9^ close the evolution system 

6 e 2 R = 2 R + 2 2 D a2 D b a { 2 - 2 A9^ , (6) 
&3«6 = 2<$ +6Wq ab , (7) 
Sl 0W = -\{6^f - <r»aW al - 8nT ab £ a £ b , (8) 

h^=^ cd qah - q \q\G lcfd H S ■ (9) 

Once initial conditions are prescribed, the only remaining in- 
formation needed to close the system are the matter term 
T ab £ a £ b in the focusing equation and q c a q d b Ci c fd£ l £ f in 
the tidal equation. Using a null tetrad (£ a , k a , m a ,m a ) 
(see Appendix [A} they can be expressed in terms of Ricci 
and Weyl scalars: 8nT ab £ a £ b _= R ab £ a £ b = 2<J> 00 and 
q c a q d b Ci c fd£ l £f = ^ortTarnh + ^ornarnt,. The complex Weyl 
scalar $q ar, d the Ricci scalar <£>oo drive the evolution of the 
geometric system ©-(O on the horizon. Being determined 
in terms of the bulk dynamics (^o relates to the near-horizon 
dynamical tidal fields and incoming gravitational radiation, 
whereas <E>oo accounts for the matter fields), fields an d $00 
act as external forces providing (modulo initial conditions) all 
the relevant dynamical information for system ©-Q on S . 

In the DH case, although the evolution system is more 
complex, the qualitative conclusions reached here remain un- 
changed. More specifically, the differential system on H gov- 
erning the evolution of 2 R is also driven by external forces 
given by a particular combination of Weyl and Ricci scalar^]. 

In the present cross-correlation approach, these dynamical 
considerations strongly support ^0 as a natural building block 
in the construction^ of the quantity hi nn (v) at 7i, to be corre- 
lated in vacuum to dPf jdu at This is hardly surprising, 
given the dual nature of a nd ^4 on inner and outer bound- 
aries, respectively. 

Particularly relevant are the following remarks. First, in 
the presence of matter, the scalar <E>oo plays a role formally 
analogous to that of ^o- Therefore, in the general case, it 
makes sense to consider <&oo on an equal footing as ^0 m the 
construction of h- mn (v). Second, Eq. © is completely driven 
by the rest of the system, without back-reacting on it. For 
this reason, although "Po (and $00) encodes the information 
determining the dynamics on the horizon, at the same time the 
evolution of 2 R is sensitive to all relevant dynamical degrees 
of freedom, providing an a ve raged response. This justifies the 
crucial role of 2 R in the construction of the effective Kf s in 
paper I. 

A serious drawback for the use of ^0 an d $00 m the con- 
struction of a quantity hi nn (v) at T-L is their dependence on 
the rescaling freedom of the null normal £ a by an arbitrary 
function on S. We address this point in the following section. 



5 In a DH, the leading term in the external driving force is indeed given by 
^o, but corrections proportional to C also appear. 

6 Constructed as in Eqs. (5J and (4) but substituting 2 R by >&o- 
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FIG. 3. Worldtubes Hi (blue) and H2 (red), respectively, associ- 
ated with two different 3 + 1 slicing {St 1 } and {£t 2 } and providing 
evolutions from a given marginally trapped surface So in an initial 
Cauchy hypersurface £rj . They illustrate the nonunique evolution of 
AHs into DHs. The foliation {(S^} (resp. {St 2 } ) by marginally 
trapped surfaces is defined by the intersections of Hi with {St x } 
(resp. H2 and |E t2 }). Note that, from the DH foliation unique- 
ness Result 1 1131 . surfaces S' t2 = Hi n E f2 are not (in general) 
marginally trapped surfaces. 

III. FUNDAMENTAL RESULTS ON DYNAMICAL 
HORIZONS 



The introduction of a preferred null tetrad on the horizon 
requires some kind of rigid structure. We argue here that DHs 
provide such a structure. We first review two fundamental 
geometric results about DHs: 

a) Result 1 (DH foliation uniqueness) II 1911 : Given a DH H, 
the foliation {S t } by marginally trapped surfaces is unique. 



b) Result 2 (DH existence) Il20l 12111 : Given a strictly stably 
outermost marginally trapped surface Sq in a Cauchy hyper- 
surface So, for eacn 3 + 1 spacetime foliation {S t } containing 
So there exists a unique DH H containing Sq and sliced by 
marginally trapped surfaces {St} such that St C S f . 

These results have the following important implications: 

i) The evolution vector h a is completely fixed on a DH (up 
to time reparametrization). By Result 1 any other evolution 
vector h' a does not transport marginally trapped surfaces into 
marginally trapped surfaces. 

ii) The evolution of an AH into a DH is nonunique. Let 
us consider an initial AH Sq C So an d two different 3 + 1 
slicings {S tl } and {S t2 }, compatible with So- From Result 
2 there exist DHs Hi = (J t St 1 and H2 = (J{ Sta> wrtn 
S tl = Hi fl S tl and S t2 — H2 H S t2 marginally trapped 
surfaces. Let us consider now the sections of Hi by {S t2 }, 
i.e., S[ 2 = Hi H S t2 , so that Hi = \J t2 S^. In the generic 
case, slicings {<S( } and {5^} of Hi are different (deform 
{S t2 } if needed). Therefore, from Result 1, S' t2 cannot be 
marginally trapped surfaces. Reasoning by contradiction, we 
then conclude that Hi and Hi are different hypersurfaces in 
Ai, as illustrated in Figf3] 



The two results above establish a fundamental link between 
DHs and the 3 + 1 approach here adopted. We denote (cf. also 
Appendix [A} the unit timelike normal to slices S t by n a and 
the spacelike (outgoing) normal to St and tangent to S t by 
s a (see Fig[T]i. We denote by N the lapse associated to the 
spacetime slicing function t, i.e., n a = —NW a t- Given a 
marginal trapped surface Sq in an initial slice So, and given 
a lapse function N, let us consider the (only) DH H given by 
Result 2. Then the unique evolution vector h a on H associated 
with Result 1 can be written up to a time-dependent reseating 
as 



h a = Nn a + bs a 



(10) 



where 6 is a function on St to be determined in terms of N and 
C [see Eq. (13[ below]. Certainly such a decomposition of an 
evolution vector compatible with a given 3+1 slicing {S t }, 
in the sense h a V ' a t = 1, is valid for any hypersurface but, in 
the case of a DH and due to Result 1, the evolution vector h a 
determined by Eq. ( fTOb has an intrinsic meaning (up to time 
reparametrization, which is irrelevant in a cross-correlation 
approach) as an object on H not requiring a 3 + 1 foliation. On 
the other hand, Eq. ([1} provides the expression of vector h a in 
terms of the null normals. More specifically, Eq. (Q]) links the 
scaling of £ a and k a to that of h a by imposing h a — > £ a as 
the DH is driven to stationarity (C -> <+ 5 t 6^ -> 0). 
Writing the null normals at H as i a = (f/2)(n a + s a ) and 
k a = (n a — s a )/ f, for some function /, expressions ([T) and 
(O for h a lead to 



pa 



N + b 



{n a + s a ) 



h a 
N 



1 



N + b 



(n a -s a ) ,(11) 



where the subindex N denotes the explicit link of H to a 3 + 1 
slicing. In order to determine b, we evaluate the norm of h a 
and note that the function C in Eq. (Q~|) is expressed in terms 
of N and b as 



C=-(b*-N 2 ) 



(12) 



On the other hand, for a given lapse N, the trapping horizon 
ShO^' = condition translates into an elliptic equation for C 
[cf. Equation ( IA12H 



- 2 AC + 2n { / )2 D c C - C 
SirT ab T a £ b 



2 



(/)(/)■* 

a ab 



(13) 



Therefore, for a given DH H associated with a 3 + 1 slic- 
ing with lapse N, Eqs. ( TOI ) and (fl2l i fix the value of b. Pre- 
scription ( fTTb provides then preferred null normals on a DH H 
compatible with the foliation defined by N. Completed with 
the complex null vector m a on St, we propose 



(£%, fcjy, m a , m a ) , 



(14) 



7 This applies, strictly, to the external part of the horizon discussed in 
Sec. Ill Al For the internal part one must reverse the evolution with respect 
to that defined by the 3 + 1 foliation: h a = -Nn a + bs a . The following 
discussion goes then through. 
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as a preferred null tetrad (up to time reparametrization) on a 
DH. To keep the notation compact, hereafter we will denote 
the preferred and simply as £ a and k a and omit the 
symbol N from all quantities evaluated in this tetrad. The 
tetrad ( TT4l > then leads to a notion of preferred Weyl (and Ricci) 
scalars on the horizon T~L. In particular, 



*o = C\ cd £ a m b £ c m d , 
*oo = ^Rab e a £ b . 



(15) 
(16) 



In summary: we have introduced preferred null normals on a 
DH H by: i) linking the normalization of l a to that of h a by 
requiring h a — > l a in stationarity; and ii) fixing the normaliza- 
tion of h a (up to a time-dependent function) by the foliation 
uniqueness result on DHs (Result 1). The latter is the rigid 
structure needed to fix a preferred null tetrad on H. In the 
particular case of constructing T-L in an initial value problem 
approach (Result 2 on DHs), the free time-dependent function 
is fixed by the lapse N of the given global foliation {£*}. 



IV. NEWS-LIKE FUNCTIONS AND BONDI-LIKE FLUXES 
ON A DYNAMICAL HORIZON 

A. News-like functions: vacuum case 

In Sec.UJwe have identified the Weyl scalar *f>o as the ob- 
ject that encodes (in vacuum and for C = 0) the relevant geo- 
metric information on the BH horizon understood as an inner 
screen. Then in Sec. |III]we have introduced a preferred scal- 
ing for on DHs. With these elements we can now introduce 
the following vectorial quantity on T-L 



with 



8tt 



dA. 



*o(v')dv', 



(17) 



(18) 



where we make use of an advanced time v parametrizing Ti 
(cf. Sec. Ill Al and Fig. [2]) and adapted to the 3 + 1 slicing at T-L 
(namely, we choose h a V a v = 2 to match the general notation 
in paper I). 

The quantity Ki could be used as a refined version of Kf 
for the correlation with dP® jdu at However, whereas 
Kf is explicitly understood as an effective quantity and, con- 
sequently, one can relax the requirement on the Af constructed 
out of 2 R in to behave mathematically as a news function, 
the situation is different for Ki in (fTTb : the geometric dual na- 
ture of ^4 and v&q would call for a newslike function character 



for A/", 



dSJ. 



Whereas expressions for the flux of Bondi momentum and 
the news function at J? + [cf. Equations (33) and (34) in paper 
I] are valid under the (strong) conditions enforced by asymp- 
totic simplicity at null infinity and in a given Bondi frame, 
no geometric structure supports the "a priori" introduction of 



quantities Ki and A^ on T-L. In particular, the news func- 
tion Af(u) is an object well-defined in terms of geometric 
quantities on sections S u C J? + , that can be expressed as 
a time integral [cf. Eq. (34) in paper I] due to the key relation 
d u N = ^4 holding for Bondi coordinate systems at J? + . On 
the contrary, the quantity A\j, defined by time integration of 
is not an object defined in terms of the geometry of a sec- 
tion S v (justifying the use of a "tilde"). Such a local-in-time 
behavior is a crucial property to be satisfied by any valid news 
function. Therefore, one would expect additional terms to "to 
(with vanishing counterparts at J? + ), contributing in J\f^ to 
build an appropriate newslike function on H. 

In the absence of a sound geometric news formalism on Ti, 

we proceed heuristically by modifying A/"4 so that it acquires 
a local-in-time character. Such a property would be guaran- 
teed if the integrand in definition ( fT8l ) could be expressed as 
a total derivative in time of some quantity defined on sections 
>S„. The scalar "Po in Eq- dTSfc does not satisfy this property. 
However, with this guideline, inspection of ^ suggests some 
of the terms to be added to [system Q-© applies to the 
EH case] so that they integrate in time to a quantity on <S„, 
namely the shear. Considering first, as an intermediate step, 
the EH case and using a tensorial rather a complex notationj, 
let us introduce a newslike tensoiQ (A^ 
ation is 



\[f )ab 



whose time vari- 



(A/f ) Qb = 



q d b c lcfd e l £ f 



that is, such that (A/"4 )o6 = -l/v^foo"^ (the global fac- 
tor 1 / V2 is required for the correct coefficient in the leading- 
order contribution). Upon time integration in Eq. ( fTSI ) and 
setting vanishing initial values at early times, this choice leads 
to 



q a b 



(19) 



1 



(20) 



If we write 



-j=j [q c a q d b C lc}d £ l t S - 
2(A4 £) ) cd (A^ ) *) C£ Vbl dv' 



(21) 



and substitute (A/"4 )ab recursively in the right hand side, we 

can express the newslike function (A/^ ) a b in terms of so 
that the lowest-order term is indeed given by expression ( fT8l . 

This identification, in the EH case, of a plausible newslike 
tensor as the shear along the evolution vector suggests the fol- 
lowing specific proposal for the newslike tensor for DHs 



1 



(/<) 



\r (n> 



(22) 



We write complex numbers as 2 X 2 traceless symmetric matrices. 
1 Note that we remove now the "tilded" notation to emphasize its newslike 
local-in-time character. 
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This proposal has a tentative character. Once we have identi- 
fied the basics, we postpone a systematic study to a forthcom- 
ing work. 



have an instantaneous meaning. However, once the newslike 
tensor M„u has been introduced in 



formal fluxes can be 
constructed by integration of the squared of these news. More 
specifically, we can introduce the formal fluxes on T~L 



B. News-like functions: matter fields 

As discussed in Sec. [II] in system ©-(O the Ricci scalar 
$00 plays a role analogous to that of From this perspec- 
tive, in the matter case, it is reasonable to define as in ( TT~8T > 



r r (l), , a m 
■A/i» = -y 



such that Ki in dT7T > is rewritten 



$00 (V) dv' , 



(23) 



dA. 
(24) 

The parameter a m is introduced to account for possible dif- 
ferent relative contributions of $o an d $oo (distinct choices 
for a m are possible, depending on the particular quantity to 
be correlated at J^ + ). However, also the function AQ, is af- 

fected by the same issues discussed above for M# , namely 
it lacks a local-in-time behavior. As in the vacuum case, we 
proceed first by looking at EHs. We then complete $oo with 
the terms in Eq. ©, so that Af^ ] (v) = -(a m /2)8 e e^. That 
is 



8nT cd £ c e d + 



so that = — (a m /2)6^\ This matter newslike function 
can be equivalently expressed in tensorial form as follows 



2V2 



'q a b ■ 



(25) 



As in vacuum, the passage from EHs to DHs is accom- 
plished by using the natural evolution vector h a along H for 
the expansion. Then, combining the tensorial form (1251 1 with 
we can write a single newslike tensor as 



Kh = — 7= 



- (a (h) 



2 



'lab 



(26) 



Interestingly, if a m = 1 the complete news tensor acquires a 
clear geometric meaning as the deformation tensor along h a , 
i.e., as the time variation of the induced metric 



1 



1 



ab 



V2 ab ' 2V2 



q a b ■ 



(27) 



C. Bondi-like fluxes on H 



dE 1 / , m) ju) 

» = ^ * AC, Af 



dv 



ab 



dA. 



(28) 



dv W 8tt 



H) 



ab 



dA. 



where their formal notation as total time derivatives is meant 
to make explicit their local-in-time nature. The purpose of 

quantities dE /dv and {dP [£,]/dv) is to provide im- 
proved quantities at H for the cross-correlation approach. In 

particular, (dP H [£,}/dv) provides a refined version of the ef- 
fective Kf s in paper I, to be correlated with (dP^ / du){u) at 
J^ + . In this context, Ki in Eq. (T7\ has played the role of an 
intermediate stage in our line of arguments. 

Of course, we can introduce formal quantities E and 
P i on H, by integrating expressions in ( fZ8b along H. How- 
ever, in the absence of a physical conservation argument or 
a geometric motivation, referring to them as (Bondi-like) en- 
ergies and momentum would be just a matter of definition^ 
Thus, we rather interpret them simply as well-defined instan- 
taneous quantities leading ultimately to a timeseries hi nn (v). 

It is illustrative to expand the squared of the news in ( l28b as 

1 r 



(H) (-H) 



ab 

~ 2 

in the 



ab 



(29) 



to be inserted in the expression for dE /dv and 

(dP [£]/dv). The relative weight of the different terms as 
we depart from equilibrium can be made explicit by express- 
ing the evolution vector as h a — £ a — Ck a [cf. Equation Q)], 



with associated crty = afJ - Ca^J and 8^ = -C9^ 



(*0 



'ab 



'ab 



ab 



[cf. Equation (1A1 1H . We can then write 



1 r 
2 



ab 



ab 



ab 



+C 2 



ab 



ab 



■(30) 



On a DH, terms proportional to a m only enter at a quadratic 
order in C. Two values of a m are of particular interest. First, 
the case a m = 0, corresponding to an analysis of pure gravi- 
tational dynamics. Second, the case a m = 1 where [cf. d27l il 



dE 



dv 



•(«) 



1 

327 



1 



{q a bfdA , 



(31) 



s v 



that admits a suggestive interpretation as a Newtonian kinetic 
energy term of the intrinsic horizon geometry. 



The motivation for introducing Kf s in paper I and Ki in 
Eq. ( TTTb [or, more generally, Ki in Eq. d24l il is the construc- 
tion of quantities on H to be correlated to quantities at J? + , 
namely the flux of Bondi linear momentum. We have been 
careful not to refer to them as to "fluxes," since they do not 



10 For instance, the leading-order contribution from matter to the BH energy 
and momentum should come from the integration of the appropriate com- 
ponent of the stress-energy tensor T a i, an element absent in j28t where 
matter contributions only enter through terms quadratic in T a i,. 
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D. Relation to quasilocal approaches to horizon momentum 
and application to recoil dynamics 

As emphasized in the previous section, the essential pur- 
pose of dE jdv and (dP [£]/dv) in (f28b is to provide 



geometrically sound proposals for hi nn (v) at T-L. Having said 
this, it is worthwhile to compare the resulting expressions, for 
specific values of a m , with DH physical fluxes derived in the 
literature. This provides an internal consistency test of the line 
of thought followed from Kf s to Eqs. (f28j. In particular, for 
a m = we obtain 



dE"" 
dv 



» 



1 

8?r 



s v 



1 

167T 



a^a^dA 



1 

167T 



s v 



WW 6 

"ab " 



dA . 
(32) 



Expression d32l allows us to draw analogies with the energy 
flux proposed in the DH geometric analysis of Refs. fil. 12211 . 
In particular, the leading term in the integrand of this expres- 
sion, <J^a^ a , is directly linked [cf. Equation (3.27) in J3l] 
to the term identified in 123l 12411 as the flux of transverse 
gravitational propagating degrees of freedoiirj. The DH en- 
ergy flux also includes a longitudinal part 1231 12411 depending 
on , absent in quantities in Eq. ( l28l i. In this sense, 

dE N jdv provides a quantity hi nn (v) accounting only for the 
transverse part of gravitational degrees of freedom I24l . l26ll27ll 
at H and therefore particularly suited for cross-correlation 
with (dPf /du)(u), which corresponds to (purely transverse) 
gravitational radiation at J? + . 

Motivated now by the resemblance of ( 1321 with the flux of a 
physical quantity, we can consider a heuristic notion of Bondi- 
like 4-momentum flux through %. Considering the (time like) 
unit normal f° to % [cf. ( lAlOl i and ( IA15H 

' -Ck a ) = -^{bn a + Ns a ) , 



(2C 



2C 



(33) 

we can introduce the component of a 4-momentum flux 
(dP^ /dv) along a generic 4-vector r/ a , as 



dP T [rf\ _ 1 



dv 



8tt 
1 

167T 



ab 



dA 



dA , (34) 



that has formally the expression of the flux of a Bondi-like 4- 
momentum. The corresponding flux of energy associated with 
an Eulerian observer n a is 

dE T dP T [n] 1 f b 

~dv~ i - V ' = dv ~ TfaJ s VT( 

(35) 



ab 



ab 



dA 



where ^= = \Jl + N 2 /2C. Analogously, the flux of linear 
momentum for £ a tangent to T, t would be 

16irJ s „ V2C 



dv 



ab 



dA. (36) 



Near equilibrium, i.e., for C — > 0, we have er^cr^) ~ C 
on DHs [cf. Equations ( 1A14M so that the integrands in expres- 
sions d35l l and (f36l > are 0(yC), therefore regular and vanish- 
ing in this limit. Considering (dP T [£\/dv)as an estimate of 
the flux of gravitational linear momentum^ through H, the 
integrated quantity P[ would provide a heuristic prescription 
for a quasilocal DH linear momentum, a sort of Bondi-like 
counterpart of the heuristic ADM-like linear momentum in- 
troduced for DHs in Ref. l28tl . by applying the ADM expres- 
sion for the linear momentum at spatial infinity i° to the DH 
section St 



P[<L\ = 



8- 



{K ab - K lab ) Cs b dA 



s, 



(37) 



In this sense, the cross-correlation methodology we propose 
here and in paper I, can be formally compared with the 
quasilocal momentum approaches in Refs. 11281 l29tl to the 
study of the recoil velocity in binary BHs mergers, showing 
the complementarity among these lines of research. 

However, attempting to derive in our context a rigorous no- 
tion of quasilocal momentum on H would require the devel- 
opment of a systematic news-functions/rameworA: on DHs, in 
particular considering the possibility of longitudinal gravita- 
tional terms as in the DH energy flux (cf. Refs. Il30l432ll for 
important insights in this topic). Such a discussion is beyond 
our present heuristic treatment, and we stick to our approach 
of considering the constructed local fluxes on T~L as quantities 
encoding information about (transverse) propagating gravita- 
tional degrees of freedom to be cross-correlated to the flux of 
Bondi momentum at J? + . 



V. LINK TO THE HORIZON VISCOUS-FLUID PICTURE 

The basic idea proposed in Ref. [01 is that certain qualita- 
tive aspects of the late-time BH recoil dynamics, and, in par- 
ticular, the antikick, can be understood in terms of the dissipa- 
tion of the anisotropic distribution of curvature on the horizon. 



11 We note that § s a^a^^dA was used in Ref. 12511 as a practical dimen- 
sionless parameter to monitor horizons approaching stationarity. Here they 
would correspond to a vanishing flow of transverse radiation. 



12 A related alternative prescription for a DH linear momentum flux would be 
given by angular integration of the appropriate components in the effective 
gravitational-radiation energy-tensor in 12411 . 
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This picture in which the BH recoils as a result of the emis- 
sion of anisotropic gravitational radiation in response to an 
anisotropic curvature distribution suggests that the interaction 
of the moving BH with its environment induces a viscous dis- 
sipation of the gravitational dynamics. The cross-correlation 
approach to near-horizon dynamics discussed in paper I and 
complemented here offers a realization of the idea proposed 
in expressing it in more geometrical terms. Indeed, the 
analysis in Sec. |IV]has led us to the identification of the shear 

and of the expansion 9^ h \ interpreted there in terms of 
newslike functions at H, as the relevant objects in tracking the 
geometry evolution. This identification permits to cast natu- 
rally the viscous-fluid picture into a more sound basis, since 
8^ and er^ have indeed an interpretation in terms of bulk 
and shear viscosities. Such dissipative features can already be 
appreciated explicitly in Eq. (0, but acquire a larger basis in 
the context of the membrane paradigm that we review below. 



A. The BH horizon viscous-fluid analogy 



Hawking and Hartle Q3M35H introduced the notion of BH 
viscosity when studying the response of the event horizon 
to external perturbations. This leads to a viscous-fluid anal- 
ogy for the treatment of the physics of the EH, fully devel- 
oped by Damour l36l l37ll and by Thorne, Price and Macdon- 
ald [38, H^], in the so-called membrane paradigm (see also 
ll40l 41]). In this approach, the physical properties of the 
BH are discussed in terms of mechanical and electromagnetic 



properties of a 2-dimensional viscous fluid. A quasilocal ver- 
sion of some of its aspects, appl ying fo r dynamical trapping 
horizons, has been developed in Il42l445ll . 

In the fluid analogy of the membrane paradigm, dissipation 
in BH dynamics is accounted for in terms of the shear and bulk 
viscosities of the fluid. The viscosity coefficients are identi- 
fied in the dissipative terms appearing in the momentum and 
energy balance equations for the 2-dimensional fluid. These 
equations are obtained from the projection of the appropriate 
components of the Einstein equations on the horizon's world- 
tube, namely evolution equations for <?W and VL a l ■ For an EH 



these equations are |4 



Miw + flW nw = 2 D a Uo + 9 —)- *d*m\ 



+8irT cd £ c q d 



(38) 



The first one [i.e., the Raychaudhuri Eq. (O not assuming a 
affine geodesic parametrization, so that ^ 0) is inter- 
preted as an energy dissipation equation. In particular, a sur- 
face energy density is identified as e = -6>W /8tt. The second 

evolution equation for the normal form Q a provides a mo- 
mentum conservation equation for the fluid, a Navier-Stokes- 
like equation (referred to as Damour-Navier-Stokes equation), 
once a momentum ir a for the 2-dimensional fluid is identified 
as 7r a = — Sl£ £ ' l /(87r) [note that f2„ is associated with a den- 
sity of angular momentum; cf. Equation dA8H . Dividing Eqs. 
by — 8it and applying these identifications we obtain 



S e e + ie >t 



87T / 107T 



Z D 



8tt 



(a {lY 



8- 



2 D 



16-k 



+ T ab tl h ■ 
- q c a T cd £ d 



(39) 
(40) 



Writing the null evolution vector as £ a — dt + V a , for 
some (velocity) vector V a tangent to St, one can write 
0« = D a V a + d t \n^q and 2o% = { 2 D a V b + 2 D b V a ) - 
O^'Qab + ^dtq a b- Then one can identify a fluid pressure 
P = «W / (8ir), a (negative) bulk viscosity coefficient £ = 



— l/(167r), a shear viscosity coefficient (j, = l/(167r), an ex- 
ternal energy production rate T ab £ a £ b and external force den- 
sity f a = —(faTcd^ ■ See also l46ll for a criticism of this 
interpretation. 

The analogue equations in dynamical trapping horizons are 
obtained from th e eq uations 5hO h and Sh^la ■ The latter can 
be written as 143144511 



S h + 6^) flW = + aVpvW* + + 2 D a { 2 D a C - 2<7ftW) + %nT ab T a h b - 6^5 h C , 

§h + 000) = 2 D a ^ - 2 D c a£ h D J^ + 8nq b a T bc r c - 6^ 2 D a C . 



(41) 
(42) 



with = ~k a h b \7 b £ a [see Eq. jA9l l. Then, by introducing a DH surface energy density e = -6» (t ) /(8tt) = 6 W /(8tt), 
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keeping the definition for ir a and introducing the heat Q a = a complete interpretation of these equations) 



1 

4w 



we can write for DHs (see 144 14511 for 



J 



^ ° {h) + ih(e (h) ) 2 + < 



rW 



T ab r a h b - 2 D a Q a - —S h C 



s h TT a + e^ h K a = - 2 D a 



8tt ) 
8tt 



16tt 



<S7T 



87, 



2 n c Uac 



8tt 



^ a 



8tt 



(43) 



We can now justify the viscosity interpretation of 0( h > and 
cr, ' by remarking that from the equations above, 9^ repre- 
sents the expansion of the fluid in the bulk viscosity term [with 
positive bulk viscosity coefficient ( = l/(167r)]. Similarly, 



a 1 ^' corresponds to the shear strain tensor and a a T b ' /{Sir) to 

the shear stress tensor. Note that a^ 1 /(8ir) and cr^ are not 
proportional in the strict dynamical case, C ^ 0, and there- 
fore one cannot define a shear viscosity coefficient /.t (in other 
words, a DH is not a Newtonian fluid). 

Finally let us consider the observer given by the (prop- 
erly normalized) timelike normal to H and let us define the 
4-momentum current density associated with this observer: 
p a = —T ao T b . Then we note that the components of p a are 
fixed by Eqs. d43l together with the trapping horizon defining 
constraint Eq. (fTTt . Indeed, p a h a = —T ao T b h a corresponds 
to the energy e balance equation, while Pbq b a = —T 0C T c q b a 
gives the momentum 7r a conservation equation, and p a r a = 
—T a bT b T a is a linear combination, using r a = 2£ a — h a , of 
the energy dissipation equation and the trapping horizon con- 
dition (5 h 8 w = 0) depending on T ab T a £ b . Given the funda- 
mental role of the latter in the geometric properties of the DH, 
in particular in the derivation of an area law under the future 
condition 8^ < 0, this suggests the possibility of using the 
component p a r a to define a balance equation for an appro- 
priate entropy density. This point echoes the discussion of a 
hydrodynamic entropy current discussed in the context of a 
fluid-gravity duality IH7l - l5lll . 



(t) 



B. A viscous "slowness parameter" 

The viscosity interpretation outlined in the previous sub- 
section allows us now to make contact with the slowness pa- 
rameter P introduced in and discussed in paper I in the 
context of BH head-on collisions. We recall that the parame- 
ter P is constructed in terms of two dynamical timescales: a 
decay timescale r and an oscillation time scale T 



P= T -. 

T 



(44) 



In our fluid analogy, the bulk viscosity term $( h > controls the 
dynamical decay, whereas the shear viscosity term cr i i s re ~ 
sponsible for the (shape) oscillations of the geometry. Given 



averaging over horizon sections we can build instantaneous 
timescales[3 at any coordinate time t as 



1 



1 



r(t) A 



(45) 
(46) 



where £| is the unit vector in the instantaneous direction of 
motion of the BH at time t. The term (£|si) in the defini- 
tions d45l l - d46b is needed for giving a timescale associated 
with a change in linear momentum [if not, we would be deal- 
ing with a timescale for a change in energy, cf. d28l)l. In other 
words, it is needed to account for the dissipation and oscil- 
lation of anisotropies in the geometry rather than for spheri- 
cally symmetric growths. This is consistent with the beating- 
frequency behavior found in the timeseries developed for the 
head-on collision of two BHs (cf. Eq. (58) in paper I). Note 
that Eqs. d43T>-(l46b provide geometric prescriptions for the 
instantaneous timescales at the merger of a binary system, 
an open problem pointed out in (H]. Combining Eqs. d44b 

and d45l-ll4"6l and denoting |crW| 2 = V 1 )^ ', we get 



P{t) 



§ St {Q Si )9WdA 
A § St {^ Sl )\aW\HA 



(47) 



As a consistency check we can verify for DHs that using 
Eqs. (IA1 II) and (IA14b and in situations close to stationarity 
(i.e., C — > 0), the following scaling holds 9^ ~ C and 
| cr C^) | 2 r*i C, so that P remains well-defined in this limit. 
For an alternative and more sound proposal for P, improving 
further the behavior when C — > 0, see Eq. JA16t . 



VI. CONCLUSIONS 

The analysis of spacetime dynamics is a very hard task in 
the absence of some rigid structure, such as symmetries or a 



their physical dimensions = [a 



ab 



[Length] " 



13 These are not the only possibility to define r and T, and therefore P, from 
viscosity scales. All variants should give though the same qualitative esti- 
mates; see Eq. )A16t . 
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preferred background geometry. However, this is the generic 
situation in the strong-field regime described by general rela- 
tivity. In this context, (complementary) effective approaches 
providing insight into the qualitative aspects of the solutions 
and suggesting avenues for their quantitative modeling are of 
much value. In this spirit, in paper I and here, we have dis- 
cussed a cross-correlation approach to near-horizon dynam- 
ics. Other interesting schemes, such as those developed at 
Caltech, that define and exploit new curvature-visualization 
tools l52ll53ll . share some aspects of this methodological ap- 
proach. 

In particular, we have argued that, in the setting of a 
3 + 1 approach to the BH spacetime construction, the foliation 
uniqueness of dynamical horizons provides a rigid structure 
that confers a preferred character to these hypersurfaces as 
probes of the BH geometry. Employed as inner screens in the 
cross-correlation approach, this DH foliation uniqueness per- 
mits to introduce the preferred normalization (fTTT > of the null 
normals to AH sections and, consequently, a preferred angular 
scaling in the Weyl scalars on these horizons. The remaining 
time reparametrization freedom (time-stretch issue) does not 
affect the adopted cross-correlation scheme, where only the 
structure of the respective sequence maxima and minima is of 
relevance in the correlation of quantities defined at outer and 
inner screens. 

Although this natural scaling of the Weyl tensors on DHs 
has an interest of its own, we have employed it here as an inter- 
mediate stage, linking the effective-curvature vector Kf (t) 
in paper I to the identification of the shear 0"*\ , associated 
with the DH evolution vector h a , as being proportional to a 
geometric DH newslike function J\f a ^ in Eq. d22l i [see also 

the role of 9^ h \ in the more general in Eq. (1261) 1. On 

the one hand, this identification provides a (refined) geometric 
flux quantity (tiPj /dv) on DH sections to be correlated to 
the flux of Bondi linear momentum (dP® /du) at (these 
DH fluxes also share features with quasilocal linear momen- 
tum treatments in the literature). On the other hand, given 
the role of and 9^ in driving the Ricci scalar 2 R along 
H [namely Eq. (f5]l and system ©-©J, the present analysis 
justifies the use of Kf s (t) in paper I as an effective local esti- 
mator at H of dynamical aspects at J' + . 

The cross-correlation analysis has also produced two im- 
portant by-products. First, we advocate the physical relevance 
of tracking the internal horizon in 3 + 1 BH evolutions. This 
follows from the consideration of the time integration of fluxes 
along the horizon and its splitting dB4b into internal horizon 
and external horizon integrals (cf. AppendixIBl. Such expres- 
sion is fixed up to an early-times integration constant, con- 
trolled by dynamics previous to the formation of the (com- 
mon) DH (and possibly vanishing in many situations of inter- 
est). Second and most importantly, from the perspective of a 
viscous-horizon analogy we have identified a dynamical de- 
cay timescale t associated with bulk viscosity and an oscilla- 
tion timescale T associated with the shear viscosity [cf. Equa- 
tions. d45l>-(l46l> and also Eqs. dA16b l- This is particularly rel- 
evant in the context of BH recoil dynamics, where the anal- 
ysis in iH shows that the qualitative features of the late-time 



recoil can be explained in terms of a generic behavior con- 
trolled by the relative values of a decay and an oscillation time 
scales. The viscous picture meets the rationale in 0] and of- 
fers an understanding of the relevant dynamical time scales 
from the (trace and traceless parts in the) evolution of the hori- 
zon intrinsic geometry, in particular, providing instantaneous 
dynamical time scales at the merger and a geometric prescrip- 
tion [cf. Equation ( l47b and also Eq. ( IA18H for the slowness 
parameter P = T/t introduced in Q§]. 

As a final remark we note that while the material presented 
here places the arguments made in 0] and in paper I on a 
much more robust geometrical basis, much of our treatment is 
still heuristic and based on intuition. More work is needed for 
the development of a fully systematic framework and this will 
be the subject of our future research. 
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Appendix A: A geometric brief 

We bring together in this Appendix the different geometric 
objects and structures that have been introduced in the text, on 
the spacetime (A-i, g a b) with Levi-Civita connection V a . 



1. Geometry of sections St 

Normal plane to S. Given a spacelike closed (compact with- 
out boundary) 2-surface S in Ai and a point p € S, the nor- 
mal plane to S, Tjj-S, can be spanned by (future-oriented) 
null vectors l a and k a (defined by the intersection between 
Tjj-S and the null cone at p). We choose a normalization 
l a k a = —1. Directions £ a and k a are uniquely determined, 
but a normalization-boost freedom remains: l' a ~ f£ a , k' a = 

Intrinsic geometry on S. The induced metric on S is given by 

q a b = fjab + k a £b + Lh ■ (Al) 

We denote the Levi-Civita connection associated with q a i, as 
2 D a . The area form on S is given by \ = ^/qdx 1 A dx 2 , 
i.e., 2 e a b = k c l d \ ca <ab, and we use the area measure notation 
dA = y/qd 2 x. 

Extrinsic geometry of S. First, given a vector v a orthogonal 
to S, we denote the derivative at S of a tensor X ai - " a " n b 1 ...b m 
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tangent to S along v a , as 

S v X ai '" an b 1 ...b m = 



(A2) 



C V X C 



di...d m ) 



where C v denote the Lie derivative along (some extension of) 
v a . Then, the deformation tensor Qw along a vector v a nor- 
mal to S 



®ab = q c a q d b^cv d = ^s v q ab , 



i 



(A3) 



encodes the deformation of the intrinsic geometry along v a . 
More generally, the second fundamental tensor is defined as 



fC c ab = q\q e b V d q c e = k c e^+e c @ ab 



('■■) 



We can express in terms of its trace and traceless parts 



(A4) 



e 



where 0W and <j^ b denote, respectively, the expansion and 
shear along v a 

ee q ab V a v b = ±S v ^q , a<$ = 0$ - q ab . 

(A6) 

Information on the extrinsic geometry of (S, q ab ) in {M., g ab ) 
is completed by the normal form Qa\ defined as 



(„) 



(A5) 



(A7) 



In particular, given an axial Killing vector <f> a on S, an angular 
momentum J[<f>\ (coinciding with the Komar angular momen- 
tum if (f> a can be extended to a Killing vector in the neighbor- 
hood of S) can be defined as 



J\<t>\ = 



1 



(A8) 



2. Evolution on the horizon H 

Given a DH H, it has a unique foliation {St} by marginally 
trapped surfaces. This fixes, up to time reparametrization, the 
evolution vector h a along H. This is characterized as being 
tangent to H and orthogonal to St, and Lie-transporting St 
onto St+st- Sht = L We write h a and a dual vector r a or- 
thogonal to T~L in terms of the null normals as 



h" 



Ck a 



(A10) 



Then h a h a = — T a r a = 2C. The expansion 000 and shear 



are written as 



(All) 



The DH is characterized by 00) = and 5 h 9^ = 0. Using 
( IA10I ) and the properties of the S v operator, the latter condition 
is expressed as 

- 2 ac* + 2n ( ^ 2 D c c - cs k e (e) = ~s e e w , (Ai2) 



an elliptic equation on C. Under the outer condition in III CI 
S k 9^ < 0, a maximum principle can be applied so that 
C > 0, with C = if and only if <^00) = (stationary case). 
Therefore, a (future outer) trapping horizon H is fully parti- 
tioned in purely stationary and purely dynamical sections. In 
other words, sections of H react as a whole, growing in size 
everywhere as soon as some energy crosses the horizon some- 
where. This nonlocal elliptic behavior is inherited from the 
defining trapping horizon condition Eq. ( IA12b . Substituting 



- 8nT ab £ a £ l 



(A13) 



into dA12l ). we recover Eq. $13[ in the text. In the spherically 
symmetric case (C = const), and using the expression for 
S e e w in (tAT3b into (IA12t . we get 



C 



ab 



(A14) 



This quantity is well-defined for any divergence-free axial 
vector cj) a . Finally, given a vector v a £ T ± S we define ll43tl 



-(») 



-k a v c S7 c £ a . 



(A9) 



Remark on S v 9^ . In <|A2} we have introduced S v in terms of 
the Lie derivative on tensorial objects. However, the evalua- 
tions of expressions such as 8 V 9^ is more delicate, since 

0« 

is not a scalar quantity on M, but rather a quasilocal object de- 
pending on S. In the general case, <5 7 „00) (with 7 a function 
on S) depends on the deformation induced on S by 7, so that 
<5 7 „0W ^ 7«5„00). This is the reason for the special notation 
S v . Properties 5 av+bw 9^ = a£ u W + b5 w 6^ (a,b G E), 
and the Leibnitz rule S v {j9^) = i6yj)9^ + ySjW still 
hold. See for instance Refs. 1I15L I20L d4I1 for a discussion of 
this derivative operator. 



3. 3+1 perspective on the horizon H 

Given a 3 + 1 foliation of spacetime {St} defined by a time 
function t, we denote the unit timelike normal to S t by n a 
and the lapse function by N, i.e., n a = — iVV a t. The induced 
metric on S t is denoted by j ab , i.e., j ab = g ab + n a n b with 
Levi-Civita connection D a . The extrinsic curvature of E t in 
M. is K ab = — j c a V c n b . We consider a horizon H, such that 
the spacetime foliation {£t} induces a foliation {St} of % 
by marginal trapped surfaces. From Result 1 in Sec. Hill this 
foliation is unique. Let us denote the normal to St tangent to 
Ef by s a . Vectors n a and s a span also the normal plane to St- 
From the condition S^t = 1 we can write h a and r a in ( lAlOb 
as 



h a = Nn a + bs a 



bn a + Ns a , (A15) 
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for some function b on St, expressed in terms of N and C in 
(lATOl as 2C = {b + N)(b - N). 



4. An improved geometric prescription for the slowness 
parameter 

In Eqs. ( 1451 ) - ( l46l i we have introduced decay and oscilla- 
tion instantaneous timescales from 0^ and , respectively, 
identified as newslike functions at H in Sec. [IV] and respon- 
sible for bulk and shear viscosities on H (cf. Sec. |Vj. This 
is not the only possibility. From the bulk and shear viscosity 
terms in Eq. (I4l1 i we define 

where K^ h ' can be expressed, in a 3 + 1 decomposition, as 

K (h) = Ns a DaN _ bs a s b Kab + Shhl f^-tl) . (A17) 



{£ a ,k a ,m a ,m a } 



Then, the slowness parameter P = Tjr in Eq. ( l44b results 

§ St {Q Si )^ h )eM)dA 



P{t) 



dA 



(A18) 



Note that, neglecting derivative and high-order terms in 
Eq. (HTT) near stationarity (C -» 0), we get k^O^CO ~ 

a ab a ^ a > so ma t f ~ 1 consistently with the expected ab- 
sence of antikick in this limit (cf. |@]). 



5. Weyl and Ricci scalars 

Let us complete null vectors £ a and k a in T^St to a 
tetrad {£ a ,k a ,(ei) a ,(e2) a }, where (ei) a are orthonormal 
vectors tangent to St- Defining the complex null vector 
m a = ^7j[( e i) a + *(e2) a ], the Weyl scalars are defined as 
the components of the Weyl tensor C a bcd in the null tetrad 



$o = C\ cd £ a m b £ c m d , 
*i = C \cd Lm b £ c k d , 
$2 = C a bcd £ a m b m c k d . 



C\ cd £ a k»m c k d , 



Ricci scalars are then defined as 



$00 = 

$n = 

$01 = 
$12 = 
$10 = 



2 R-ab^ & i 

-jR ab (£ a k b + m a m b ) 

-\R ab £ a m\ 

^R ab k a m b , 

-lR ab £ a m b , 



$4 = c\ cd m a k b m c k d , 



(A19) 

$21 = -±Rabk a m b , 
$02 = -7>R ab m a m b , 

$22 =-\Rabk a k b , 

$20 = - \Rabrn a m b , 
K=±R. 

(A20) 



Appendix B: Relevance of the 3+1 inner common horizon 

In this appendix we emphasize the role of the inner hori- 
zon present in 3 + 1 slicings of BH spacetimes, and discussed 
in Sec. Ill Al when considering the time integration of fluxes 
along the DH history. This is of specific relevance to the dis- 
cussion made in Sec. II V Dl but it also applies to more general 
contexts. 

Given a flux density -^g- (17, v) through % of a physical 
quantity Q(v), we can write 



Q( v ) = Q( Vo ) + 

= QW+ / F c (v')dv' , 



(Cl,v')dA ) dv' 
(Bl) 



whereH F c [v) = § Sv sign(C)^j(w')dA This requires a 
good parametrization of H by the (advanced) coordinate v, as 
well as an initial value Q(vq). Finding such an initial value 
is in general nontrivial and this is precisely the motivation to 
consider in this section the evaluation of the fluxes along the 
whole spacetime history of H, though from a 3+1 perspective. 

Given the 3 + 1 slicing {S t }, we can split the integration 
along the DH into an external and an internal horizon parts, 
as discussed in Sec. Ill Al Denoting by v c the advanced time 
associated with the moment t c of first 3+1 appearance of the 
horizon, T~L is separated into the inner horizon "Hint labeled by 
Vq < v < v c and the outer horizon H ou t labeled by v c < v < 

^:H = H in t U H oxt =([j V0 < v < Vc 5.) U (U c <„<oo Sv 
We can then rewrite Eq. (IB 1 b as 



Q(v) = Q(v )+ F c {v')dv' = Q(v )+ I F^(v')dv'+I F^(v')dv' 

J Vq 

p2v c -v 

= Q(v )+ F^(2v c -v")dv" + I F^iv'W. 



(B2) 
(B3) 



The sign sign(C), +1 for spacelike and —1 for timelike sectors of H, 
corrects the possibility of integrating twice (null) fluxes through "H, when 
timelike parts occur in the world-tube of the trapping horizon H. Note 



that sign(C) appears under the integral since a section S v can be partially 
timelike and partially spacelike, i.e. the evolution vector h a can be timelike 
or spacelike in differemt parts of S v . 
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where Fq* and denote, respectively, the flux of Q along 
the internal and external horizons. Note that in the second 
term in ( IB3b we have inverted the integration limits in order 
to have an expression which is ready to be translated for an 
integration in t. 

The coordinate v is not usually adopted in standard 3 + 
1 numerical constructions of spacetimes. Because of this, 
we employ the time t defining the slicing {St}. Although 
the t function is not a good parameter on the whole H, it 

I 



correctly parametrizes the evolution of both the inner Hint 
and outer H cx t horizons separately: H = Hint U H e *t = 
(Ut>t c s t nt ) U (Ui>t S t xt ) ■ Considering the splitting in 
Eq. (IB2b . the use of t in the flux integration is perfectly valid 
as long as the t-integration includes both the standard external 
horizon part and an internal horizon part. 

From Eq. (IB3b we write 



Q(t) = Qo + J F^(t')dt' + J F^\t')dt' 

= Qo+ f F^{t')dt' + f F^\t')dt' + Res(t) , 



(B4) 



where Qq is a constant and the error Rcs(t) 

/oo 
F£\t')dt', (B5) 

must be taken into account, since we cannot integrate up to 
t — >• oo during the 3 + 1 evolution. This error satisfies 



Res(t) — y as t — y oo, so that the evaluation of Q(t) by 
ignoring Res(i) in Eq. (IB4b improves as we advance in time 
t (cf. Figure|2]). Of course, this approach requires a good nu- 
merical tracking of the inner horizon, something potentially 
challenging from a numerical point of view (see lf55ll for a 
related discussion). 
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